Nernst effect of Dirac fermions in graphene under weak magnetic field 
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The derivation for the transport coefficients of an electron system in the presence of temperature 
gradient and the electric and magnetic fields are presented. The Nernst conductivity and the trans- 
verse thermoelectric power of the Dirac fermions in graphene under charged impurity scatterings 
and weak magnetic field are calculated on basis of the self-consistent Born approximation. The 
result is compared with so far the available experimental data. 
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I. INTRODUCTION 

The Nernst effect as well as the thermoelectric power 
is a sensitive probe of the impurity scatterings in an elec- 
tron system. Recently, the Hall and the Nernst effects in 
graphene have been studied experimentally at relatively 
strongir— and moderate^ magnetic fields. Graphene in 
most of the experiment devices is absorbed on the sur- 
face of Si02. There are strong evidences that the charged 
impurities in the substrate are responsible for the carrier 
density dependences of the electric conductivity^"— and 
the Hall coefficientii as measured in the experiments by 
Novoselov et al.t^ At strong magnetic field, the carriers 
are in the Landau quantized states. In the interior of the 
system under the strong magnetic field, the carriers are 
mostly localized around the charged impurities. Though, 
the Hall effect seems weakly dependent of the impurity 
scatterings because the current is most likely conducted 
by the edge states that are not localized The standard 
Green's function theory of many body system has diffi- 
culty to treat the charge transport under scatterings of 
the charged impurities in a strong magnetic field since it 
deals with the bulk states of the electrons. On the other 
hand, at weak magnetic field when the effect of Landau 
quantization is negligible, the standard Green's function 
theory should be applicable for investigating the magne- 
tothermoelectric transports of the electron system. 

Based on the self-consistent Born approximation 
(SCBA) for Dirac fermions under the charged impu- 
rity scatterings, we have recently developed an elec- 
tronic transport formalism for graphene j^^ii^ It has 
been shown that the experimentally measured electric 
conductivity, the inverse Hall coefficient^^ and the ther- 
moelectric power—"— are successfully explained by our ap- 
proach. In this work, along the same approach, we study 
the Nernst effect of Dirac fermions as a function of the 
carrier density under a weak magnetic field. Though 
there exists no experimental measurements of the Nernst 
effect in a weak field so far, we show that our obtained re- 
sults for the transverse thermoelectric power could qual- 
itatively compare with the experimental measurements^ 
in magnetic fields of moderate strength. We intend to 
examine to what extend the theory is valid in dealing 
with the transport properties of graphene. 



Meanwhile in doing this work, we present a derivation 
of the transport coefficients of an electron system under 
the temperature gradient and the electric and magnetic 
fields being applied. 

The model of electrons in graphene is established from 
its energy band structure in the first Brillouin zone 
corresponding to a honeycomb lattice. At low car- 
rier concentration, the low energy excitations of elec- 
trons in graphene can be viewed as massless Dirac 
fermions. ^^'^-"^- That is, the energy linearly depends on 
the momentum around the two Dirac points in the first 
Brillouin zone. Using the Pauli matrices cr's and r's to 
coordinate the electrons in the two sublattices (a and h) 
of the honeycomb lattice and two valleys (around the two 
Dirac points 1 and 2), respectively, and suppressing the 
spin indices for briefness, the Hamiltonian of the system 
is given by 



H = '^Tpl[vk ■ ar^ 



kq 



>fe (1) 



where = (cLi 



"fca2) fermion oper- 



ator, the momentum k is measured from the center of 
each valley, v (~ 5.86 eVA) is the velocity of electrons, 
/i is the chemical potential, V is the volume of sys- 
tem, and Vi{q) = ni{—q)vo{q) is the charged impurity 
potential. ^^'^^ Here, ni{—q) is the impurity density and 
VQ{q) is given by the Thomas- Fermi (TF) type 



(2) 



where qxF = '^kpe^/ve is the TF wavenumber, kp = 
y/Trn (with n as the carrier density) is the Fermi 
wavenumber, and e ^ 3 is the effective dielectric con- 
stant. For briefness, we hereafter use units oi v = h = 
^^(the Boltzmann constant) = 1. 

With SCBA^i22 ^j^g Green's function 

G{k,uj) = [uj + ^ — k ■ (jTz ~ 'S{k,(jj)]~^ 
= go{k,uj) + gcik,uj)k ■ aTz 

and the self-energy Y,{k,ijj) of the single particles are 
determined by coupled integral equations.— The diago- 
nal and off diagonal parts, go and gc respectively, of the 
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Green's function can be expressed as 

5o,c(fc,w) = [sr+(fc,cj) ±.g_(fc,a;)]/2 

with g± as the upper and lower band Green's func- 
tions. Corresponding to the SCBA to the self-energy, 
the current vertex correction Tx{k, wi, a;2) is given by the 
ladder-diagrams approximation as shown in Fig. 1(a). 
Tx{k,uJi,uj2) is expanded as 

3 

{k,UJl,UJ2) = ^y,(fc,a;i,W2)AJ(fc) (3) 

3=0 

where A^(k) = r^cr^, Af{k) = a^a ■ k, A^{k) = a ■ ka^, 

A^{k) ~ ■ kuxd ■ k, and yj{k, 011,012) are determined 
by four-coupled integral equations<^ The functions yj de- 
scribe how the current vertex is renormalized by the im- 
purity scatterings from the bare one Aq(A:). 

II. FORMALISM 
A. General formula of the transport coefficients 

To study the Nernst effect of graphene, we consider 
the electronic transport of Dirac fermions under weak in- 
plane temperature gradient VT, electric potential 0, and 
weak magnetic filed B = V x A perpendicular to the 
graphene plane. Here, A is the vector potential. Since 
the temperature gradient VT is not a dynamic quan- 
tity, we cannot directly apply the linear response theory 
(LRT) to treat the current response to VT. To use LRT, 
one usually introduces a fictitious gravitational potential 
(that couples with the Hamiltonian) following the work 
of Luttinger and obtains the transport coefficients using 
the Einstein argument relating the currents response to 
the external perturbations.^'^ Here, we present a deriva- 
tion from a microscopic point of view. 

First, following the idea of Luttinger^ suppose the 
system with a variable temperature T(r) in locally equi- 
librium everywhere in space. Specifically, consider the 
system is divided into small cells but microscopically 
large enough. The Hamiltonian of the cell at rj with 
chemical potential fJ-{rj) and temperature T{rj) are given 
by Hj = f^^^^_^drip''{r)[h{r) - /x(rj)]-0(r) with h{r) = 

(TTz • (— iV + A) + Vi{r) and tli{r) as the operator of Dirac 
fermions in real space. The distribution function of the 
system is then given by 

pi = Z-'exp[-J2H,/T{r,)] (4) 
j 

where Z is the normalization constant. Instead of con- 
sidering T(r), we intend to find out an equivalent system 
determined by an effective Hamiltonian H^f / at constant 
temperature Tq. Its distribution function is 

P2 = Z-iexp(-i/e///To). (5) 



From pi = p2, we have H^ff = ^^j?o/T(rj). Sup- 
pose each cell is macroscopically so small enough that 
the summation can be replaced with integral over space. 
Heff reads 

Heff = J dfii\r)[h{r) ~ ^(r)] o T^i'ir) 

= J dfip^r)[h{r) + h{r) o vl'(r) + <I>(r) - Ho]ip{r) 

with h{r) o ^{r) = {/i(r), *(r)}/2 (here {A,B} is the 
anti-commutation relation between A and B) and 

$(r) = po-Wr)/T{r), (6) 
*(r) = To/T(r)-l. (7) 

Here, Tq and /iq are the average temperature and chem- 
ical potential, respectively. In the limit T(r) — >■ Tq, we 
have [^(r), $(r), 'l'(r)] — > [^o,0, 0]. By so doing, the sys- 
tem with variable temperature T(r) in the local equilib- 
rium state is now described by an equivalent one under 
the potentials [^{r), '^{r)] at constant To and fiQ. 

Now we go back to the original problem: How do the 
currents respond to the temperature gradient VT? Ini- 
tially the system is in the equilibrium state of -ff[0, 0]. 
With gradually turning on VT, the system H[0, 0] be- 
comes unstable because the equilibrium state shifts to 
i?[$(r), 4'(r)] and thereby currents are produced. {In 
the shifting process [0,0] [$(r), ^'(r)], Tq and ^0 are 
kept as constants.} The perturbation here is the differ- 
ence -ff[0,0] - i7[$(r), *(r)]. This is different from the 
usual case that the perturbations are due to applying the 
external dynamic potentials and the initial equilibrium 
state given by H[0, 0] keeps unchanged. 

Generally, in addition to VT, with the external scalar 
potential </> being applied, the system under consider- 
ation is 0]. With respect to the equilibrium sys- 
tem the perturbation is 0] - iJ[$(r), *(r)]. 
Mathematically, we have 

H[(j),0] = H[^,^+ H[(l),0]- H[^,^ 
= H[<i>,^]+H' 

with H' given by 

H' = j dfiii\r)[-h{r) o ^'(r) + <j){r) - $(r)]V'(r) 

= J dfV'^(r)[$i(r) + ^(r) o $2(t')]V'(0 

where ^i{r) = (j){r) - $(r) - ^Q'i!{r), ^2{r) = ~*(''), 
and ^(r) = h{r) — jiQ. Here $i(r) and ^2{t) take the role 
as the perturbation potentials. In the limit VT 0, the 
negtive forces V<l'i(r) and V^i(r) read 

V<I>i(r) = \I[(j){r)- ^i{r)]=eE, (8) 
V$2(r) = -ToV[l/T(r)]. (9) 
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Hereafter we denote fiQ and Tq simply as /i and T, re- 
spectively for briefness. According to LRT, we need to 
find out the corresponding currents determined from the 
equations of continuity. We here consider the relevant 
currents. 

(i) For the potentials $i and $2 (actually the corre- 
sponding vector potentials), the coupling currents to be 
determined are Ji(r) and Jiir), respectively. One may 
consider to convert the facts tp^{r)ip{r) (coupled to $1 in 
H') and tp^ {r)£^{r)'tlj{r) (coupled to $2) into the respec- 
tive currents in the picture of if [$, \E'] (from which the 
perturbed system evolve). The resulted currents then 
contain the terms of $1 and $2. However, since H' is 
already linear in $1 and $2, we just only need to con- 
sider them in the picture of H[0, 0]. For the unperturbed 
system H[0,0], they are the particle current and the heat 
current. 



J2{r) = J^{r) ^tl;\r)i{r)ojtl){r). 



(10) 
(11) 



with j ~ Stz- 

(ii) The currents of the reference (using subscript r) 
system 5'] itself can be obtained from the known 
results of the particle current Jir and energy current 

in Ref. [l^. The results are 

Jir{r) = J{r)o[l + -^{r)], 
J2r{r) = J^{r)- iiJir{r) 

= J^{r) o [1 + 2*(r)] + [$(0 + ^^'(r)] o J(r). 

Their averages under vanish. 

(iii) The currents of the system iJ[(^,0] under physical 
(using subscript p) observation are 

Jivir) = J>), (12) 
J2p{r) = j'^{r) + cp{r)oJ{r). (13) 

In terms of Jir{r) and J2r{i'), and $1 and $2, they read 

Jip(r) = Jir(0 + 'A(r)o$2(r), 

J2p{r) = J2,(r) + Ji(r)o$i(r)-H 2/2(0 o $2(0. 

The observed current densities are their averages in the 
system 0]. 

By the linear response theory, we have 

(Jip(r)) = (Ji(r)o$2(r-))o-i^"-V$i-i^i2.v$2^ 

(14) 

{J2p{r)) = (Ji(r)o$i(r))o+2(J2(r)o $2(0)0 

. v$i • V$2, (15) 

where (• • • )o is the average under -ff [$, ^], K'^^^s are 2x2 
constant tensors [with respect to the directions of the co- 
ordinates hereafter denoted by subscripts (a;, y) or (/i, v), 



the superscripts 1 corresponding to particle current and 
2 to heat current] determined by the Kubo formula 



IC^ = - lim Imx'^(r2 + iO)/17 



(16) 



with x*-' (51 -|- iO) as the retarded response function. In 
the bosonic Matsubara frequency fi™, x*"* reads 



1 

v 



dT^^^^^^TrMT)Jj{0))o (17) 



where Ji — J drJi{r). Now that {Jip{r)) given by Eqs. 
(ITil) and are already explicitly linear in the pertur- 
bations $1 and $2, in calculating the averages for the 
constants in Eqs. (|14p and (|15p . the equilibrium point 
can be shifted to [0,0]. To the first order of V$i, 
the term -(/(r) o <I>j(r))o's in Eqs. JTll) and are 
calculated as 

(Jl(r) o $^(r)>o = (Jl(r) o (f - V$j))o. 

For the diagonal elements {Jix{r) o x)o = M*^, we have 

which can be shown by expanding V'('') in terms of 
the eigen states of ^(r). For the off-diagonal elements, 
(Axir) o y)o = AT^y, we get 



m: 



xy - {Jix{r)y)o, 

= {[J^xir)y-J,yir)x]) 0/2, 



(18) 



because the system is invariant under rotation around 
the z-axis perpendicular to the plane. The element M'^y 
is the magnetization of the electrons. Substituting the 
results into Eqs. (|14l) and dTSl) . we get 



(Jipir)) 
{•hpir)) 



with 



-TV" • V$i - TV^^ . (19) 
JV^i • V$i - Ar22 . V$2, (20) 



(21) 
(22) 
(23) 
(24) 













- M 






- M 


^22 


- i^22 


- 2M^ 



These forms have been obtained in Ref. |2J with a phe- 
nomenological approach using the Einstein argument and 
the consideration classifying the transport components in 
each kind of the local currents as the observable currents. 

Though we consider the Dirac fermions here, the 
derivation above is valid for general electron systems. 

For a non-interacting system such as the one consid- 
ered here, the average (• • • )o under -ff[0,0] reduces, in 



principle, to the independent single-particle problem. Us- 
ing the bases of single particle states {\n)} for a given 
impurity configuration, one can obtain formally the ex- 
pressions for N^^ and N^"^ = N'^^ as 



(a) 



TV" 

^12 



du) 
2^ 



(25) 



— f{^)u:C{u:), (26) 

ZTT 



with 



= ^ejnmjninS{l^ - " Cm))i (27) 

n^rn 

where f'{uj) — df{uj)/duj with f{uj) the Fermi distribu- 
tion function, jnm — is the eigenvalue of ^(r) 
of the state \n) and (• • • )i in Eq. ^7} means the average 
over the impurity configurations. These forms have been 
obtained in Ref. [2^. For the readers' convenience, we 
give a derivation in Appendix. The compact form given 
by Eqs. ([^5)) is so obtained because of large cancella- 
tions between K^^ = K^^ and M. The tensors N^-^ and 
N^'^ = N'^^ are thus related via the function (^(a;). 

However, Eq. (j27p is not a convenient formula to start 
with. To proceed, one needs to derive C{uj) in terms of 
the Green's function. From Eq. (|17l) . carrying out the 
r-integral, we have 

Xllmm) = '^''J rf^"Tr(G(r,r',iw„)Ji^ 



xG{r',r, iujn + 

— Pfiiy {cOn ; + ^m) 



(28) 



where G{r,r' ,iujn) is the Green's function for a given 
impurity distribution, a;„ is the fermionic Matsubara fre- 
quency, {■■ ■)i again the average over the impurity con- 
figurations, and the function P^y(u;i,a;2) is so defined 
by the equation. Taking the analytical continuation 
Q,rn + iO, we have for iV", 

/duj 
^f'iu;)Re[PiLj-,uj+) - P{lo- .uj-)] 



where uj^ — ui HO, and 



dz-^[P{z',z)^P{z,z%,=, 



Therefore we get 

C{uj) = Re[P(w",w+) -P(w",w") +f(a;+)] 







FIG. 1: (Color online) (a) Self-consistent Born approximation 
for the Current vertex correction, (b) Diagrams for calculat- 
ing the function Pxy{ijJ-i,u)2)- The vertex j is associated with 
the vector potential Aj. fc* = fc ± q/2. 



B. Hall and Nernst conductivities of Dirac 
fermions 



The Nernst effect describes the response of the trans- 
verse current to a temperature gradient VT in the pres- 
ence of a perpendicular magnetic field B but i? = 0. It is 
reflected by the Nernst conductivity iV^^. Here, we study 
it in the limit of S 0. 

In the limit of _B 0, the elements A^^^ and N^.'^ are 
independent of the magnetic field. They are related to 
the electric conductivity a = e^iV^^ and thermoelectric 
power S = —N^'^/eTN^^l which we have given in previ- 
ous works j -^^i-^^ Though the element N^y was calculated 
previously for studying the Hall coefficient, the function 
Cxy{uj) was not given explicitly.^^ To calculate N^y, we 
here need to find out the explicit expression for Cxy{uj). 

As in the calculation of the Hall conductivity in the 
limit of _B — Ofii by introducing the vector potential via 
A{r) ~ A{q) exp{iq-f) with B = iqx A{q) and taking the 
limit q ^ 0, one obtains in terms of the average 

of the multiplication of three current operators i^^i^'' As 
shown in Fig. 1(b), Pxy{uJi,uj2) is obtained as 



A; 

x[G{k+,LU2)Ty{k+,uj2,u:i)V{k+,k~,LJi) 

+ V{k+ ,k- ,LJ2)ry{k- ,UJ2,UJl)G{k- ,LJl)]} 

where the factor 2 stems from the spin degeneracy, 
/c± = k±q/2, V{k+,k-,uj) = G{k+ ,uj)f{k+ ,k- ,u,uj) ■ 

AG{k^ ,Lu), and Ta:{k,uji,uj2) = Tx{k,k,LL>i,uj2) as given 
by Eq. ([3]). The vertex rp(fc~ , fc+, wi, ^2) satisfies the 
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4x4 matrix equation 

ki 

xr^{kY,k+,uji,uj2)G{k+,uj2). (29) 

To find out the limit of g ^ 0, we need to expand tlie 
right hand side of Eq. (|28|) to the first order in q and 
then use B = iq x A. The manipulation is tedious but 
elementary. We only outline the key points in the deriva- 
tion below. 

(i) The expansions of the Green function G{k^,uj) 
and the vertex functions r^(fc^, wi, a;2) can be easily ob- 
tained by definition. The most involved expansion is 
for the vertex function Tx{k^ ,k^ ,u!i,u!2)- By expand- 
ing both sides of Eq. ([29|) to the first order in q, one 
gets rf,{k^,k+,uji,uj2) = Tf,{k,uJi,uj2)+it,(k,uJi,uj2)-q/2 
with 7^(fc, wi, 012) determined by 

7p(/c, wi,a;2) = ^njWQ(fc - fc')G(fc', wi)7^(fc',a;i,a;2) 



k' 

X [VG(A:', wi)r^(fc', wi, W2)G(fc', c^a) 
-G(fc',wi)r^(fc',wi,W2)VG(fc',W2)] 



(30) 



where V means the gradient with respect to k' . 
(ii) From the identity 



-^^7i,[(Vfc + VfcOvo(fc- fc')]Tr[G(fc,wi)7^(fc,wi,W2) 



kk' 



xG{k,u2)G{k' ,u2)T,{k' ,uj2,uji)G{k' = 0, (31) 

performing the integral by part in the left hand side 
of Eq. (PT|) and using Eq. (15(11) and the equation for 
T^{k,uji,uj2), we obtain 

^^'^{'lt,(k,uji,i^2)[^kG{k,uj2)T^{k,U2,uji)G{k,uji) 

k 

+G(A:,cj2)r^(fc,W2,wi)VfcG(fc,wi)]} 
= - ^ Tr{[VG(fc, c^i)r^(A:, c^i, t^2)G(fc, UJ2) 

k 

-G{k, uji)T^{k, wi, a;2)VG(fc, uj2)]VkTu{k, uj2,uji)}. 

(iii) For V{k^,k^,uj), using Eq. 1^^, we get the ex- 
pansion 

G{k+,u)r^{k+,k-,uj,uj)G{k-,u) 



— — G(fc, a;) -|- iazA ■ \b(k, uj)k(j) — a(k, uj)(f)k] ■ q. 
oka 



(32) 



where (p is the angle of fc, is the unit vec- 
tor in (f) direction, and the first term in the right 



hand of Eq. (1321) comes from the Ward identity 
G{k,uj)ra{k,uj,uj)G{k,uj) = dG{k,uj) / dk^- The coeffi- 
cients a{k,iij) and b{k,uj) are determined by solving Eq. 
([50)) . Since the final result depends on their combina- 
tion a{k,uj) + h{k,uj) = c{k,uj), the function c{k,uj) = 
z{k,uj) + [gQ{k,uj) — g'^{k,uj)]X{k,uj) is determined by 
following equations: 

z{k,uj) = [gQ{k,oj)gc{k,u!) - go{k,u!)g'^{k,uj)] 
x[yQ{k,cj,oj) - y3{k,uj,uj)] 
-gc{k, uj){go{k, uj)[yQ{k, u, uj) + ysik, cj, cj)] 
+'2gc{k,u!)yi{k,u!,u!)}/k, 

X{k,Lo) = l^n,«2(fc~fc'){z(fc',w) 

A;' 

+g+{k',Lo)g^{k',uj)X{k',uj)}, 

with g'{k,uj) = dg{k,uj)/dk. 

Using the results given above, one gets a final expres- 
sion Pxy{u}i,u}2) = [Z{u}i,L02) - Z{uj2,uji)]/2i with 

2Be^ ,dG2 ar,,2i 

^(0.1,0.2) = ^LTr[Gir..i2(^^^ 



V 

dG2 dry^2i . 



ic{k, uJi)Tx^i2G2Ty^2icrz 



where Tx,ij — Tx{k,uJi,LLjj) (and the same meaning for 
Ty^ij), Gj — G{k,LL!j). Since Pxy{uj^ ,uj^) = 0, we have 
Cxy{u)) = ImZ(w^,a;+) -|- YieYxy{io^) , and 

dzlTiiR{z) 

-00 

with R{uj) = g§7[^(w',w) - Z(w,w')]c.'=c.- 

Knowing the function Gxyii^)-, we can calculate N^y 
and the Nernst conductivity N^y according to Eqs. ([25]) 



and (|26p , respectively. Since we are interested in the low 
temperature cases, we here give their expressions in the 
limit of r ^ 0. They are 



TV. 



Nt 



(33) 



1 f^d 
— ImZ(0",0+)- / -^Rei?(iz) 

^Im[|-Z(c.-,c.+)U.o + i?(0+)], (34) 
6 ouj 



In obtaining Nj^y, the use of the expanding Cxy{<^) ~ 
Cxyi'^) + tjjC'^y{Q) has been made. In addition, the z- 
integral in Y{.eYxy{()) reduces to the path- integral around 
the negative axis (—00,0). Because R{z) is an analytical 
function (by definition) in the upper and lower z-plane, 
respectively, the integral path has been deformed to the 
imaginary axis, giving rise to the last term in Eq. p3p . 
The integral along the imaginary axis is simple to handle 
for the numerical calculation since there is no singularity 
in the Green's function. If this term is neglected, the 
expression for the Hall conductivity axy = e^TV^^ will 
reduce to the form as in the previous work?li The fact 



that the contribution from this term is neghgible wih be 
checked later. In analogous to the thermoelectric power 
S = E^/{VT)^, we define S^y = E^/{VT)y which de- 
scribes the production of the transverse electric field due 
to a temperature gradient in the absence of current fiow. 
By setting {Jip{r)) — in Eq. ([T9|) . Sxy is obtained as 



S:,y = -{Saxy + eN^l/T)/a. 



(35) 



Clearly, S^y includes two parts. The quantity Saxy de- 
scribe the process for the response of the transverse elec- 
tric field to VT through the way: Because there exist 
a longitudinal electric field (with the magnitude propor- 
tional to the thermoelectric power S) due to the tem- 
perature gradient, the current could fiow transversely by 
the Hall process. In different to this, eN^y/T implies 
an additional transverse field due to a transverse current 
response directly to VT. 
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III. RESULTS 

The functions gQ^c{k,uj) and yj{k,uji,uj2) and their 
derivatives with respect to k and lo are involved in the 
calculation. In a recent work,— we have described how 
to numerically solve the corresponding integral equations 
to determine these functions. In our numerical calcula- 
tion, we take the impurity density as rii = 1.15 x 10~'^a~^ 
(with a as the lattice constant of graphene) the same as 
in our previous works 

The numerical results for the Hall conductivity axy and 
the inverse Hall coefficient = Ba^ juxy as functions 
of the carrier concentration 8 (doped carrier per carbon 
atom) are shown in Fig. 2. The calculations with and 
without the last term in Eq. ([55]) for both results of R^^ 
and Oxy (normalized by gB in inset of Fig. 2) are almost 
indistinguishable. This fact means that the contribution 
from the last term in Eq. (1551) is negligible small. The ex- 
perimental data^^ (symbols) and the classical prediction 
= —nec are also plotted in Fig. 2 for comparison. 
As we have stated previously, the divergence of R^^ at 
6 — stems from the vanishing of dxy while the conduc- 
tivity a remains finite. 

At low T, the Nernst conductivity N^y is proportional 
to T^. In Fig. 3, we depict eNiyx'^B as a function of the 
electron doping concentration 5. Because of the electron- 
hole symmetry, N^y is even for 5 — > —5. N^y comes 
mainly from the first term in the square bracket in Eq. 
(p4| . This is similar to the case as in N^y . The derivation 
dZ{uj^ ,uj^) /duj\i^=o is very delicate. Here, it cannot be 
considered simply as cx daxyip) / because the scatter- 
ing potential here depends strongly on the electron dop- 
ing concentration. For comparison, we also depict the 
result for Suxy/TB in Fig. 3. In the regime of 5 studied 
here, the magnitudes of both quantities are overall the 
same. But at low 5, eN^l/T^B is biger than Suxy/TB. 
In the absence of the electric field, the transverse current 
density is solely determined by Jx = —N^y{S/T)y/T. 



FIG. 2: (Color online) The inverse Hall coefficient (in 
unit of 10~"^Tesla/Ohm) as a function of S. The solid (red) 
and dashed (blue) lines correspond respectively to the calcu- 
lations with and without the contribution from the last term 
in Eq. (|33|) [Inset, the corresponding results for axy normal- 
ized by aB with B is in unit of Tesla (T)]. Dot-dashed line: 
classical result. Symbols: experimental data (Ref. il2, ). 



In Fig. 4, we show the transverse thermoelectric power 
Sxy divided by TB as functions of 6 at the average impu- 
rity densities rii = 1.15xl0~^a~^ andrii = 1.6x lO'^^a"^. 
Sxy is a superposition of Saxy and eN^y/T. As a result, 
Sxy is linear in T and B. At large carrier doping, both 
of Saxy and eN^y/T have about the same contribution 
to Sxy While at low doping, eN^y/T is predominant. 
Sxy changes sign at low doping because eN^y/T does. 
The factor Saxy + eN^y/T of Sxy decreases quickly as 
d decreasing at low S regime. Its slop normalized by a 
negative constant (~ —a/S, a oc 6 at large 5 and is flat at 
very low doping) is almost the behavior of Sxy The dip 
in Sxy corresponding to the maximum of the slop. So far 
there exist no measurements of the Nernst effect of Dirac 
fermions in a weak magnetic field. But for the purpose to 
make a qualitative comparison, the experimental results 
for the transverse thermoelectric power of Ref. 2 with the 
magnetic fields B — 1 Tesla (T) (squares) and 3 T (dia- 
monds) are plotted. Here the ratio of the first non-zero 
Landau level/the Fermi energy is about 0.35 for a doping 
level atl.OxlO-''andS = lT. Though the strengths 
of these magnetic fields could not be regarded as weak, 
the data indicate a tendency toward to the theoretical 
prediction as the magnetic field is decreasing. 

The present calculation is based on the assumption 
that the impurities are randomly distributed. Actu- 
ally, at low carrier doping, graphene is an inhomoge- 
neous system due to the impurity correlations in the 
substrate as observed by experiment i^"— It seems there 
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FIG. 3: (Color online) The coefficient eNll/T^B (red line 
with circles)(in unit of fj,VK~^T~^e^ /h with h the Planck 
constant) as a function of S. The dashed line (green with 
squares) is Sa^y/TB for comparison. 

exist the electron and hole puddles. Nonetheless, in 
each puddle, the average number of the carrier is less 
than unity. Moreover, the mean free path of the elec- 
trons is much longer than the length scale (~ tens 
nanometers) of the puddles. Within a mean free path, 
a carrier can transfer through many such puddles. The 
puddles can be thus regarded as the microscopic wrin- 
kles. In addition, at low carrier concentration, there ex- 
ists significant quantum coherence between the upper- 
and lower-band states resulting in the minimum elec- 
tric conductivity,—!^"— the unconventional behaviors of 
the inverse Hall conductivityiiii^ and the thermoelectric 
power Therefore, the carrier must be treated quan- 
tum mechanically and the present approach seems plau- 
sible. 

Finally, we compare our calculation with the semi- 
classical Boltzmann theory. By the Boltzmann theory 
within the relaxation time-r approximation, the function 
gk describing the difference between the disturbed distri- 
bution function and the Fermi distribution function / is 
determined by^ 

9k = T^v-{%VT + eE) + Te{vxB)-Vk9k 
oik T 

(36) 

using units of c = ?i = 1 again. Here — vk — for 
electrons in graphene. From Eq. (1361) . one obtains 

Nil = T'rn/S. 
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FIG. 4: (Color online) The coefficient Sxy/TB of transverse 
thermoelectric power (in unit of fiV K~^T~^) as a function 
of 5. The solid (red with small symbols) and dashed (green) 
lines are the calculations at impurity densities iii — 1.15 x 
10~^a~^ and 1.6 x 10~^a~^, respectively. The squares (blue 
at B = 1 Tesla) and diamonds (cyan at B = 3 Tesla) are the 
experimental results of Ref. 2. 

The inverse Hall coefficient is Ba'^/axy — ne. The 
Matt relatior*^ is given by = T'^-n'^a' /ie^ with 

cr' = tc^/tt. However, using the Mott relation, one ob- 
tains zero Nernst conductivity Nly — because axy is 
constant independent of the chemical potential fx. This 
is different from the present result that the Nernst con- 
ductivity is in the order of Sa^yT / e as shown in Fig. 
3. 



IV. SUMMARY 

In summary, we have derived the formula for the trans- 
port coefficients for the Dirac fermions in graphene in the 
presence of the temperature gradient, the electric field 
and the magnetic field. The derivation is valid for gen- 
eral electron systems. It is different from the usual per- 
turbation process with the external dynamic potentials 
applied (in that case the original equilibrium state is un- 
changed) since the perturbation due to turning on the 
temperature gradient shifts the equilibrium state. The 
physical observed system is in the perturbed state with 
respect to this equilibrium state. 

On the basis of self-consistent Born approximation, we 
have studied the Nernst effect of the Dirac fermions un- 
der the charged impurity scatterings and weak magnetic 
field in graphene. The transverse thermoelectric power 
is closely related with the Hall conductivity and the lon- 
gitudinal thermoelectric power for which the theory has 
been shown to be in good agreement with the experi- 
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ment. The Nernst conductivity is dealt with the similar 
approach as for the Hall conductivity. The present cal- 
culation is a prediction to the Nernst conductivity of the 
Dirac fermions in graphene under a weak magnetic field. 
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APPENDIX 

We here give a derivation of Eqs. ([25)) and (1251) . Con- 
sider the response function xab{t) defined by 



1 



X^^(r) = --(T.A(r)i3(0))o 



(37) 



have 



1 



-7rRe(a„mfem„)(5(fi + - im)\li 
-TrRe{anmbmn)5{fl + C« - S.m)])i- (40) 

where P means taking the principle value in the sum- 
mation, and in the P term the exchange n ^ m and the 
use of Im(a„„6„„) = -Imial^bl^^) = -Im(a„„6„„) have 
been made in the last equality. Substituting the result of 
Eq. (gni) into Eq. we obtain 



K 



AB 1 / FT (- u \ W" fmjP 

= T7^Z^ [Im(a„„6„„)— — — 2- 
-7rRe(a„m6m„)/,'j(5(^„ - Cm)])i 



(41) 



where A and B (of Hermitians) can be either the current 
or the heat current operator. Using the basis of the single 
particle states we have 



Air) = J dfe^[°''']>^(r)a(r)V'(r)e-^[°'"]" 



where a„m = {n\a\m) and cj^ (c„) creates (annhilates) a 
particle in state \n). The function xab{t) is written as 



,AB 



V 

n 

where /„ = (cJjC„) is the Fermi distribution function. 
Taking the Fourier transform, we get 



AB 



^ ^ I ^ ^ 0^nmbmn}i- (39) 



where /3 = 1/T. Taking the analytical continuation 
i^i — > 17 -I- iO [for which the first term in the braces 
in Eq. (|5^ can be disregarded because of fi^ > 0], we 



where = dfn/d^n- 

Note first, the desired forms given by Eqs. ([25]) and 
([25]) are obtained from the contribution from the last 
term in the square brackets in Eq. (HTj) . 

Second, for the diagonal elements, the first term in 
the square brackets in Eq. (|4T|) vanishes. In the present 
case, A and B are vectors. For the diagonal elements. 



cx bx 



and cix,nmb. 



is real. 



Therefore, in following, we will consider only the con- 
tribution from the first term in the square brackets in Eq. 
(PT|) for the case of off-diagonal elements. Denote it as 



(/n frii)P 
V ^ {in - U)' 



Im(a„m6m„) (42) 



dropping the symbol of the average (• • • ); for briefness. 
We only need to prove that the off-diagonal element of 
R{A, B) is cancelled by the corresponding matrix element 
ofM'^y given by Eq. 

(i) A = B = J. Using J = i[$(r),rl, R{J, J) reads 

R{JxiJy) — ^ ^ {fn /m)Im(^nmymn) 

n 

= 0. 

(ii) A = J and B ^ J'^ with jQ{r) = {j(r), ^(r)}/2. 
Using = j„m(Cn+Cm)/2 = ir„„ (^^ _ ^2^)/2, one gets 

R{Jx, Jy) = ^ E + - /m)Im(a;„„ ?;„„). 
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On the other hand, we note last equahty. Therefore, we have 

Mxy — ^ ' fn[jx,nmymn jy,nmXmn] Ii{Jxj M^y = — ^ ^ (/wC'i frnCm)^^^i'^nmymn) 

nm nni 



nm 



— 21^ ^ ^ ^^^^ ^ /m)(^n CrM)Ini(^nmymn)5 ^ 

(hi) The case of A = and _B = J is the same as (ii) 
where we have made use of the exchange n ^ m in the and N^'^ — iV^^. 
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